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We present a dynamical description of slow relaxation processes based on the extension of On- 
sager's fluctuation theory to systems in local quasi-equilibrium. A non-Markovian Fokker-Planck 
equation for the conditional probability density is derived and from it, we obtain the relaxation 
equation for the moments. We show that the fluctuation-dissipation theorem can be formulated in 
terms of the temperature of the system at local quasi-equilibrium which is related to that of the 
bath by means of a scaling factor revealing lack of thermal equilibrium. Our theory may be applied 
to a wide variety of systems undergoing slow relaxation. We discuss in particular slow dynamics in 
glassy systems and Brownian motion in a granular gas. 

I. INTRODUCTION 

Systems whose global evolution results from the existence of a wide variety of time and length scales frequently 
exhibit slow relaxation dynamics, which is manifested through the aging behaviour of the correlations and the non- 
exponential decay of the response function Q-Q. To characterize the dynamics of these systems is one of the most 
challenging problems of nonequilibrium statistical mechanics nowadays. Many efforts have been made in this sense 
by proposing different approaches including spin models, stochastic dynamics and thermodynamical approaches (see, 



for example, Refs. |4j|-|l3l|'). 

Experiments performed in diverse systems as in amorphous polymers 14] and supercooled liquids and glasses Q], 



lixpcr 

Q-E3, 



, seem to indicate that these systems undergo, in general, non-Markovian and non-stationary dynamics leading 
to the peculiar behaviour of the two-time correlation functions and to the dependence of the relaxation times on the 



. 

initial time 



In this paper we present a dynamical description of slow relaxation systems based on the formulation of a Fokker- 
Planck equation for non-Markovian processes. This equation is obtained through the generalization of the Onsager 



o 
o 
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theory, originally proposed to account for the dynamics of fluctuations around equilibrium states |18| and extended 

• to nonequilibrium stationary states |l9|. to nonequilibrium aging states. 

_ _ i 

An important feature of the dynamics of these systems is that at the time scales of experiments, they still relax 
and never reach thermal equilibrium with the bath. We have assumed the existence of a local quasi-equilibrium state 
|2f| characterized by a non-stationary probability density which can be used to determine average quantities as, for 
example, the local quasi-equilibrium temperature Jjjl ■ The entropy of the system has been expressed in terms of the 
conditional probability density by means of the Gibbs entropy postulate, and used to derive the generalized Fokker- 
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Planck equation by taking into account probability conservation and the rules of nonequilibrium thermodynamics 
[2^ 1 . For stationary processes, this equation was derived in j^] and j^, and has subsequently been treated in the 
literature of stochastic processes in (23- 13- 

The unified description of slow relaxation systems we propose can be applied to, among others, glassy systems 
granular matter and anomalous diffusion problems Q • 

The paper is distributed as follows. In Sec. II, we will analyze the non-stationary dynamics of slow relaxation 
systems under the framework of a generalized Onsager theory. Section III is devoted to the formulation of relaxation 
equations and of the generalized fluctuation-dissipation theorem for non-stationary systems. An illustration of the 
theory presented is given in Sec. IV by analyzing the dynamics of systems having two well-separated time scales 
(glassy systems). In Sec. V, we will discuss the Brownian motion in a granular gas. We will compare our results with 
those obtained by means of kinetic theory and simulations. In the discussion section we will summarize our main 
results. 



II. FOKKER-PLANCK EQUATION FOR NON-STATIONARY DYNAMICS 

The non-stationary dynamics of slow relaxation systems can be characterized at the mesoscopic level by means of 
the conditional probability density P(a , to\a, t), which satisfies the normalization condition J P{a 0l to\a,t)da = 1. 
This quantity depends in general on a set of fluctuating variables a — (ai, ...,a n ), or fields a(fc), |3.1 | . necessary to 
univocally determine the state of the system during the relaxation and which may represent, for instance, fluctuating 
positions and velocities, order parameters, reaction coordinates or hydrodynamic fields. 

The evolution in time of P{a , to\a, t), with a the state of the system when the fluctuation takes place at to (in 
glassy systems literature it is referred to as waiting time or time elapsed after quenching), is governed by the continuity 
equation 

d d r 1 

— P(ao,io|a,t) = -— • |P(so>*o|",%«(a!o,a;*o,£)J , (1) 

expressing probability conservation during its natural evolution in a-space. The probability current 
P(&o-> to\cx, t)v a (a , a; to, t), with v a {a , a; to, t) a stream velocity in a-space, has to be determined in order to render 
the dynamical description complete. 

This task can be accomplished by assuming that the scheme of nonequilibrium thermodynamics [2] remains valid in 
a-space. Thus, a linear relationship can be established between the diffusion current P(io|t)w Q (tQ, t) and its conjugate 
force in a-space, which can be obtained from the entropy production calculated in the framework of Onsagcr's theory 
[2^. To derive the entropy production we will assume that, for sufficiently long times t, the system reaches a local 
quasi-equilibrium state, characterized by the entropy per unit mass s qe (t) and the probability density P qe (o., t), related 
formally to the conditional probability through lim^oo P(a , to\a, t) = P qe (a, t). 

In accordance with the variation of the entropy s qe (t) satisfies the Gibbs equation 



T(t)5s qe (t) = 5e qe (t) - J fi qe (t)5P(t)da, (2) 

where the temperature T(t) of the system at local quasi-equilibrium will be assumed to be only a function of time, 
e qe (t) is the mean internal energy, and (J. qe (t) the chemical potential at local quasi-equilibrium, which do not depend 
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on the a's. For simplicity we have omitted the explicit dependence of the probability density on phase variables. The 
entropy production a corresponding to the diffusion process in a-space can be computed by using the Gibbs entropy 
postulate |2J 



*(*) = -— / P(to\t)ln^lda + Sqe (t) : (3) 

m J Pqe{t) 

where s(t) is the entropy of the system, ks the Boltzmann's constant and to the molecular mass. Notice that, according 
with mesoscopic nonequilibrium thermodynamics (MNET), the entropy s qe (t) appearing in Eq. I|3|l. constitutes a 
macroscopic quantity whose variations are governed by Eq. • The evolution in time of s qe (t) can also be described 
with the MNET scheme by using the expression of Eq. J3J), and the corresponding Eq. @ referred to the equilibrium 
state (sec, for example, 23| and in that case, by taking the variation 5s qe — s qe ({P qe + 5P}) — s qe ({P qe }) at 

first order in SP and combining the resulting expression with Eq. we obtain 



P qe (t) = P e e*B m TW^ e <*\ (4) 

which incorporates the topology of the energy landscape through e qe , and where P e is the equilibrium probability 
distribution. 

By combining the rate of change of s qe (t) obtained from Eq. @ with the time derivative of Eq. @, after using 
Eq. and integrating by parts assuming that -P(iol0^a(^Oji) vanishes at the boundary, we obtain |34| 



dt 



s(t) 



1 d 



T{t) dt 



e(t) 



T(t) 



P(t \t)v a {t Q ;t) 



k B T(t) P(t \t) ; 



da, 



(5) 



where Se(t) = 5e qe {t) 



k B T(t) 



J P(tg\t)S hi P qe (t)da is the variation of the nonequilibrium internal energy e(t). If no 



external thermodynamic forces as, for example, velocity or temperature gradients act on the system, it may be shown 
by using the Gibbs-Duhcm relation that 5e(t) — 0. As a consequence of this, the entropy production is only due to 
the underlying diffusion process in a-space. Moreover, in the last term of Eq. 10) one may identify the nonequilibrium 
"chemical potential" /j,(to;t) as 



= + (6) 

TO Pqe{t) 



Important to notice is the fact that )i(to]t) may depend, in general, on interaction potentials or entropic barriers 
through the corresponding dependence of P qe (t). 

According to the rules of nonequilibrium thermodynamics expressing currents as linear functions of the conjugate 
thermodynamic forces 22] , from the entropy production JSJ) the velocity w Q (to ; t) can now be expressed in terms of 
obtaining 



£,.(/„:/! = -BUu:l)- \X(^t) + ^^-^\nP(t \t) 

to oa 



where we have identified the generalized force 



0, Q 



(7) 
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(S) 



m 

and introduced the coefficients B_(t ;t) which may, in general, depend on the state variables a. These coefficients, 
incorporating memory effects through their time dependence, can be determined in terms of the conditional relaxation 
functions of the corresponding variables a, j^]- Notice that, in general, the -By's do not satisfy reciprocal relations 
since for non-stationary systems and non-conservative interactions, the time reversal symmetry and the principle of 
detailed balance are not necessarily fulfilled. The time dependence of the B^ 's is the origin of anomalous diffusion, 
which is characterized by the power law behaviour in time of the correlation function j^J] , \<>\ , and the nonexponential 
time decay of the relaxation functions [J. 

By substituting the expression of the probability current obtained from Eq. into Eq. Q, we obtain the 
Fokker-Planck equation 



Y( APU \+\ i k B T(t) dP(t \t) \ 

X(a,t)P(t \t) + — ( 9 ) 



dt da 

Similar equations describing stationary non-Markovian dynamics were obtained in Refs. 26]-30], to study linear and 
nonlinear transport. From Eq. ©, it follows that the characteristic relaxation times of the system (to;t), are 
functions of the waiting time to as a consequence of the non-Markovian nature of the dynamics. 

At the time scales considered, the system relaxes towards the time-dependent local quasi-equilibrium state and 
then thermal equilibrium between system and bath can not be established. The complete description of the dynamics 
can then be performed by giving the relation between the system and bath temperatures. In order to derive it, the 
influence of the bath is incorporated through the second law Sst = Ss qe + Ssb > 0, with st the total entropy of the 
system and sb that of the bath. From this expression, by assuming that the variation Sst takes place at constant 
volume (see Eq. {SJ), it can be expressed as Sst = [1 — c] j^jj[Se qe (t) — 6g qe (t)], where we have taken into account the 
Gibbs equation (J2J, used the Gibbs free energy per mass unit Sg qe (t) = J fi qe (t)SP(t)da and defined c = — ^f- > 1, 
with ob the entropy production of the bath. The equality occurs when the process is reversible, i.e., ob = —o- In a 
similar way, by using the Gibbs equation of the bath we obtain Sst = [1 — c 1 ] j^jy<5es(i), with es{t) and T B (t) the 
energy and temperature of the bath. A simple comparison of the expressions for Sst leads to the relation 



T(t) = AT B (t) , (10) 



where we have defined A 



,_(t)-8e q e(t) 



Since 5g qe (t) = s qe 8T{t) — fJ, qe Sp, we conclude that the factor A also 
depends on the density and the temperature of the system. When the process is reversible, A = 1 since thermal 
equilibrium between system and bath is attained. In this case, P qe (t) reduces to a local equilibrium probability 
density in which the temperature is that of the bath 



35], 36]. 



III. RELAXATION EQUATIONS AND THE FLUCTUATION- DISSIPATION THEOREM 



Once the Fokker-Planck equation has been derived, we will proceed to establish the relaxation equations for the 
moments of the probability density. The analysis of the general case can be accomplished by considering that the 
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coefficients are functions of the state variables [3^ . For the sake of simplicity, we will focus on the case in which the 
coefficients Bij are only functions of time. 

First, we will consider the two-time correlation function 

C ai a 0l {t Q ,t) = J a t aojP(ao'*o)-P(ao'*o|a, t)da da, (11) 

whose evolution equation may be found by taking the time derivative of Eq. I|ll|) , substituting Eq. Q and integrating 
by parts, obtaining 



d_ 

dt 



C ai u j(to,t) — —Bu(to;t) (Xi(t)a j(to)) , 



(12) 



where we have defined (Xi(t)a>oj(to)) = / Xi(a,t)aojP(a ,to)P(a ,to\a,t)da Q da. This equation implicitly provides 
an operational definition of the transport coefficients Bij, showing aging effects through their dependence on the 
waiting time. The relation between the dissipation of the energy and the response of the system to an external time- 
dependent force e(t) — eQ(t — to), with Q(x) being the step function, may be determined by means of the evolution 
equation for the first moment of the probability density, (ai(t)) — J ctiP(a 0l to)P{a , to\a, t)da da. In the presence 
of the external force, the system is described by means of the corresponding Fokker-Planck equation 



9P(to\t) _ d 

— a*— - to= (<0,t) 



X^t)P(t \ t) + ^^A dP ^ 



da 



~-k(t)P(t \t)], 



(13) 



where the external force has introduced the additional drift e(t) . According to linear response theory [^3 j the response 
Ra iaj {to, t) of the system is then given through the functional derivative 



(*o,t) 



6(ai{t)) 
Sej(t ) ' 



(14) 



Notice that, due to the time dependence of the coefficients S«(tg; t), the response function will contain modifications 
to the Markovian response. 

The response function can be calculated from the equation for the first moment (cti{t)) which follows from the 
Fokker-Planck equation l|13fl . We obtain 



R ai a 0j (to,t) - ^^y^ C ' Q l Q oj(*o,i), (15) 

which constitutes the formulation of the fluctuation-dissipation theorem. This expression indicates that the 
fluctuation-dissipation theorem remains valid when the system is in local quasi-equilibrium; the temperature en- 
tering that expression is that of the system at local quasi-equilibrium. If one uses Eqs. (|10fl and (|15l) . it can also be 
expressed as 

m A^ 1 d 

R aiOt0j (t ,t) = — — —C^fat). (16) 



The result (|16|) can in general be interpreted as a violation of the fluctuation-dissipation theorem Q and |j| . When 
analyzing its equivalent form (| 1 5J1 it expresses validity of the fluctuation-dissipation theorem in terms of the local 
quasi-equilibrium temperature. 

The quantity A has been derived in Section II from thermodynamic grounds. An expression in terms of correlation 
functions can be obtained from the equal-time correlation function $> a a(t) = \ /a • aP(a Q ,to)P(a ,to\a,t)da Q da. 
After substituting Eq. 10 and integrating by parts, one obtains the evolution equation 

^ aa (t) = -B(t ;t) (X(t)-a(t)) + ^I^B(t Q ;t), (17) 
at m 

where for simplicity we have considered B_(to;t) = B(to;t)l, with 1 the unit tensor. For times t 3> l3^ 1 (to;t), the 

time derivative of Eq. I|17l) may be neglected and the system relaxes to the reference state, leading to the expression 



AS MiM&® (18) 
where {■■■) qe represents the average in respect to the local quasi-equilibrium probability density P qe (t). This equation 
indicates that A may in general be expressed by means of the moments of the probability density. According to Eqs. 
(|10|l and (|18|l . it follows that the temperature T(t) entering the Gibbs equation J2J results from the contribution of 
the energies related to the corresponding degrees of freedom a by 



T(t) = j^(X(t).a(t)) qe . (19) 

An important fact to be noticed from this expression is that the temperature is not robust since it implicitly depends 
on the number of relevant variables or observables used to describe the system, as has been evidenced in |38j. Given 
that the probability density of the reference state depends on the initial condition, it follows that T(t) is also a 
function of to. In the case when the force is linear, Eq. (|19|l reduces to T(t) = ^ (a 2 (£)) ?e , which coincides with the 
corresponding result obtained when local thermodynamic equilibrium between system and bath exists. In that case, 
when a = u with u the velocity of a particle, Eq. I|19[) becomes the kinetic definition of the temperature. 



IV. SLOW DYNAMICS IN GLASSY SYSTEMS 



As a first application of the theory presented, we will analyze the dynamics of systems characterized by two sets 
of variables: fast (3- variables and slow a- variables. The (3- variables relax to a local equilibrium state characterized 
by the temperature of the bath Tb, whereas the a- variables relax to a quasi-equilibrium state characterized by the 
temperature T(t). 

The probability density is now a function of the state vector T = (a, (3). Since [3- variables equilibrate with the 
bath, we may assume that the probability density characterizing the reference state can be factorized in the form 
Pqe{i) — Pie{l3)P qe {(^, t) , with Pi e {f3) and P qe {o_, t) stationary and non-stationary distributions. Consequently, s qe (t) = 
s ie + s ge where Tg Ss^ and T(t) Ss qe enter Gibbs equations similar to Eq. Consistent with the analysis of Sec. 
II, the probability densities of the reference states are given by 

P le (f3) = e ^-^- e '<=(«] and p qe (a,t) = P e e^w l ^~ eqA - )] . (20) 
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By taking the time derivative of Eq. © and using Eq. (|20|l and the corresponding continuity equation, the entropy 
production can be written as 



P(t \t)v{T,t) ■ -^Q(t ;t)dT, 



(21) 



where we have defined the nonequilibrium Massieu function 



efo ; t) s ^lnPfo|t) + f^> , ^ 
m T B T(t) 



(a) 



(22) 



which is the appropriate function in the entropy representation 



H, 0. 



Since v(T, t) 



and 



or 



_a_ _o 

df3 ' dc 



— ) , from the entropy production i|21|) we can obtain the expressions for 



v_p and v a (see Eq. |7|). In general, the rules of nonequilibrium thermodynamics introduce couplings between the 
velocities and v a and the corresponding forces through the mobility B_ which is in general a function of ioj t, QL 
and f3. 

In this case, since no couplings appear due to the separation of time scales, the mobility matrix B_ can be written 
as the direct sum of two submatrices 5 a and Bf, For simplicity, we will assume that the B^'s are constant since the 
/3-variables decay to a local equilibrium state |22j. The Bf^s will be functions of two times, B % = Bfj(to,t). 

Now, by taking into account that (3- variables thermalize at Tb and a- variables at T(t), after substituting 

the expressions for and w a into the continuity equation one obtains the generalized Fokker-Planck equation 



dP(t \t) 
dt 



oa — 



X a (a)P(t \t) + 



k B T{t) dP(t \t) 



m 



<:)a 



+ 



d_ 

df3 
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Xf{(3)P{t \t) 



k B T B dP(t \t) 



d(3 



(23) 



where, according with Eqs. © and (|2*U|) . Xf(f3) = 



de lc ((3) 



andr(a) 



de qs (a) 



00 — V— I ~ da ■ 

One may obtain the Fokker-Planck equation for the slow a- variables by multiplying Eq. (|23|l by P(T_Q,to) and 
integrating over and f3 taking into account that the current P(to\t)v_p vanishes at the boundaries. Our result is 
Eq. @ . In similar form, the evolution equation for the (3- variables can be derived by integrating over a a and a. The 
resulting Fokker-Planck equation is 



dP%t) 
dt 



op = 



k B T B d 



(24) 



m d(3 — 

where P(/3,t) is the single-time distribution function. 

The computation of the moments of the probability density following the procedure indicated in Sec. Ill, allows 
one to find the statistical expression of the temperatures 



K B \ — —I le 



(25) 



and 
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T(t) = £-(X«(a,t)-a) qe . 



(26) 



Eq. I|25|) implies that the fast degrees of freedom satisfy the equipartition law and its energy is a quadratic function in 
/?'s, whereas (|26|1 is a definition of the temperature of the system in terms of the correlations at the quasi-equilibrium 
state. 

^From Eq. (|23[) . we may compute the response of the system to an external perturbation and derive the fluctuation- 
dissipation theorem related to both sets of variables. Following the procedure indicated in Sec. Ill, the response 
functions of the fast variables satisfy the relation 



TTi 1 d 
KB -LB Oto 



(27) 



where Rp i p aj {ti3,t) = ^§^y- The corresponding expression for the response of the a-variables is given in Eq. ||TSJ. 
From Eq. I|27|) it follows that the fast degrees of freedom satisfy the fluctuation-dissipation theorem in terms of the 
temperature of the bath. 



V. BROWNIAN MOTION IN A GRANULAR GAS 

As a second application of our formalism, we will discuss in this section the Brownian motion of a grain in a granular 
gas, which has been previously addressed in by means of kinetic theory for systems with inelastic collisions. The 
phase-space vector is now a = (m, r) , where u and r are the velocity and position of the Brownian grain. The 
corresponding continuity equation for P(to\t) is given by 



dP(t \t) 
dt 



V • OP(to\t) = • [P{t a \t)Mto]t)] 
ou 



(28) 



We will assume that the bath is in a non-stationary homogeneous cooling state characterized by the kinetic tem- 
perature Ts(t), jslj. The kinetic temperature is given by Tit) = jr^ f \P(uo,to)P(uo>to\u,t)duodu, and has to be 
understood as the temperature of a dilute gas of Brownian particles suspended in the granular gas. According to Eq. 
(0J, in this case the local quasi-equilibrium probability density can be taken as 



P qe (t) = e 



(29) 



where /j, qe normalizes P qe {t). 

Once the P qe (t) has been specified, the Fokker-Planck equation for the Brownian gas may be derived by using Eqs. 
©, 1128(1 . 129|l . Following the steps indicated in Sec. II, we finally obtain the generalized Fokker-Planck equation 



dP(to\t) ^ d 
-^t— + V - uP= d-u B ^ t) 



uP(t \t) H A — — — 



(30) 
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where we have assumed locality in it-space for which the velocity vs(to; t) couples only to the gradient of the chemical 
potential at the same value of u in accordance with Eq. ©. By considering the action of an external force field 
e(i) which can be an impulsive one (tapping force), we may derive by following the steps indicated in Sec. Ill, the 
fluctuation-dissipation relation 

m A^ 1 d 

Ruu (t ,t) = ———C m (to,t). (31) 

Recent computer simulations of diffusion in dilute granular gases confirm the validity of Eq. I|31|l , [4 1| . The corre- 
sponding integrated version of the generalized fluctuation-dissipation theorem has been also obtained in dense granular 
systems in Now, by comparing Eq. 13Ufl with the corresponding one found in jjyl by following kinetic theory 
methods, we may identify 

A= (l + A) . BM = i(l + A)7e(*o,*), (32) 

where the restitution coefficient A(< 1) of the Brownian grains has been assumed constant, and 7 e (io> t) is the friction 
coefficient for the case when elastic collisions take place, |3l| . 

^From the Fokker-Planck equation i|30|) , it follows that the evolution equation of the second moment of the velocity 
of the Brownian grains can be expressed in the form 

jT(t) = -2B(t Q ; t) [T(t) - AT B {t)] . (33) 

which is a particular case of Eq. I|17|l . Eq. i|33[l establishes that the relation between system and bath temperatures 
at the local quasi-equilibrium state given by Eq. (|10|l . is valid for times t 3> -B _1 (io;^)- 

The relation between T(t) and T B (t) can be generalized in order to describe shorter time scales by incorporating 
the thermal relaxation of the heat bath through the evolution equation j^T B (t) = — C(t)T B (t), with £(t) the cooling 
rate. Assuming that < ((t) < 2B(tQ\t), the time evolution of is then given by [31| 

/V) +2(l-OB(t ;t)P^ = 2B(t ;t)A, (34) 



dtT B {t) K s; v u ' >T B (t) 

2B(t :t) 



where we have defined £ = 2 B(t\) • From Eq. (|34H . we then obtain 



T W=(T^) r ^W' ( 35 ) 

which gives the relation between the system and bath temperatures for times t 3> [(1 — £) B(to; t)] . 

By contracting the description, we enter the diffusion regime described by the density p(f > ,to\r,t) = mj P(to\t)du, 
whose evolution is governed by a Smoluchowski equation. This equation may be derived by approximating the 
hierarchy of evolution equations of the first three moments of the probability density P(to\t), [43j. One obtains, 

J^>(*o|t) = D eff {t a ,t)V 2 p{t \t), (36) 
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where D e ff(to,t) is the effective diffusion coefficient given by 



^ . ( 37 ) 



_ k B T B {t) A 

mB(t ,t) (l-£) 

which has been obtained by taking into account Eqs. (|34|) and l|35l) . Eq. (|37(l was derived in jsjj by means of kinetic 
theory. In particular, when the cooling rate of the bath satisfies ((t) <C 2B, the expression of the diffusion coefficient 
(|37|l reduces to D e ff(to,t) = m ^^ t ) -A- In the general case, the dynamics of the system may be described by means 
of the corresponding Fokker-Planck equation in which the probability density of the reference state has the form 
Pqe{t) — e k B T B(t) A ^ £)I^« e 2 U 1^ anc [ ^- ne coefficients ^4 and B(to,t) satisfy more general relations than Eq. (|32ll . 

The Fokker-Planck and Smoluchowski equations derived previously may also describe diffusion in a colloidal system. 
In that case, the diffusion coefficient ()37|) accounts for the violation of the Stokes-Einstein relation reported in dense 
granular systems and supercooled colloidal fluids 



VI. DISCUSSION 

In this paper, we have proposed an extension of Onsager's theory, originally formulated to describe nonequilibrium 
aged systems, in order to analyze slow relaxation dynamics. Our theory is based on the assumption of the existence of 
a local quasi-equilibrium state in phase space characterized by a non-stationary probability density. This hypothesis 
entails the existence of a time-dependent temperature of the system T(t) reflecting lack of thermal equilibrium between 
system and bath, which enables us to formulate a Gibbs equation. The relation between the local quasi-equilibrium 
and bath temperatures was given in terms of a scale function obtained by taking into account the influence of the 
bath. 

The entropy production derived from the Gibbs entropy postulate leads to the explicit expression of the probabil- 
ity current in phase space and, consequently, to the generalized Fokker-Planck equation @. Memory effects were 
incorporated in the description through the time dependence of the coefficients entering into that equation. They 
are responsible for the dependence of the characteristic relaxation times on the initial (waiting) time. Our formalism 
could straightforwardly be extended to the case of a description in terms of fields a(k), [3], Q]. 

We have shown that the fluctuation-dissipation relation remains valid in terms of the local quasi-equilibrium temper- 
ature T(t) of the system whose value depends on the time scale considered. The fact that the local quasi-equilibrium 
temperature depends on time and on the number of variables relevant to characterize the state, clarifies why two 
systems with slow relaxation dynamics put in thermal contact do not necessarily satisfy a "zero law" . 

As a first application, we have analyzed the slow dynamics of glasses. We have shown that the dynamics of the 
system can be described by a generalized Fokker-Planck equation governing the time evolution of both the fast and 
slow degrees of freedom. By using this equation, we have verified that the fluctuation-dissipation theorem (|27|l is 
satisfied by the fast modes. For the slow modes the theorem <|15[1 is valid in terms of the local quasi-equilibrium 
temperature. 

We have also studied the problem of Brownian diffusion in a dilute granular gas reproducing the Fokker-Planck and 
Smoluchowski equations derived previously by means of kinetic theory of granular gases. These equations may also be 



11 



used to describe the dynamics of supercooled colloidal suspensions leading to a generalization of the Stokes-Einstein 
relation. 

The theory presented offers a unified description of systems undergoing slow relaxation and constitutes a general- 
ization of the Onsager theory describing the evolution of fluctuations around noncquilibrium aging states. 
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